Biophysical Journal - 

Biophysical Society 

Article 


A Bidomain Model for Lens Microcirculation 

Yi Zhu , 1 Shixin Xu, ■ Robert S. Eisenberg, ' and Huaxiong Huang ■ 

’Department of Mathematics and Statistics, York University, Toronto, Ontario, Canada; 2 Centre for Quantitative Analysis and Modelling, Fields 
Institute for Research in Mathematical Sciences, Toronto, Ontario, Canada; department of Applied Mathematics, Illinois Institute of 
Technology, Chicago, Illinois; and department of Physiology and Biophysics, Rush University, Chicago, Illinois 

ABSTRACT There exists a large body of research on the lens of the mammalian eye over the past several decades. The objec¬ 
tive of this work is to provide a link between the most recent computational models and some of the pioneering work in the 1970s 
and 80s. We introduce a general nonelectroneutral model to study the microcirculation in the lens of the eye. It describes the 
steady-state relationships among ion fluxes, between water flow and electric field inside cells, and in the narrow extracellular 
spaces between cells in the lens. Using asymptotic analysis, we derive a simplified model based on physiological data and 
compare our results with those in the literature. We show that our simplified model can be reduced further to the first-generation 
models, whereas our full model is consistent with the most recent computational models. In addition, our simplified model cap¬ 
tures in its equations the main features of the full computational models. Our results serve as a useful link intermediate between 
the computational models and the first-generation analytical models. Simplified models of this sort may be particularly helpful as 
the roles of similar osmotic pumps of microcirculation are examined in other tissues with narrow extracellular spaces, such as 
cardiac and skeletal muscle, liver, kidney, epithelia in general, and the narrow extracellular spaces of the central nervous sys¬ 
tem, the “brain.” Simplified models may reveal the general functional plan of these systems before full computational models 
become feasible and specific. 

INTRODUCTION 

Biological systems require continual inputs of mass and en¬ 
ergy to stay alive. They are open systems that require the 
flow of matter and specific chemicals across their bound¬ 
aries. Unicellular organisms can provide that flow by diffu¬ 
sion to and across cell membranes. Diffusion is not adequate 
over distances larger than a few cell diameters, i.e., larger 
than 2 x 10 -6 m, to pick a number. For that reason, multi¬ 
cellular organisms cannot provide those flows to their cells 
by diffusion itself. Multicellular organisms depend on con¬ 
vection to bring materials close enough to cells so diffusion 
to and across cell membranes can provide what the cell 
needs to live. 

The circulatory system of blood vessels—arteries, 
veins, and capillaries—provides the convection in almost 
all tissues. But there is one clear exception: the lens of 
the (mammalian) eye. The lens does not have blood ves¬ 
sels, presumably because even capillaries would seri¬ 
ously interfere with transparency. The lens is large, 
much larger than the length scale on which diffusion it- 
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self is efficient. The lens must provide nutrients through 
another kind of convection, a microcirculation of water 
that moves nutrients into the lens and rinses wastes out 
of it. The microcirculation is in fact driven by the lens 
itself, without an external “pump.” The lens is itself an 
osmotic pump. 

The lens is an asymmetrical electrical syncytium in which 
all cells are electrically coupled one to another, with a nar¬ 
row extracellular space between the cells (see Fig. 1). The 
extracellular space is filled with ionic solution in “free diffu¬ 
sion” with the plasma outside cells. It may also contain 
specialized, more or less immobile proteins and specialized 
polysaccharides, as well as containing obstructions formed 
by the connexin proteins themselves. The intracellular space 
behaves very much as a large single cell would, with the bio¬ 
ions of classical electrophysiology (Na + , K + , Cl - ) free to 
move without much resistance from cell to cell and many 
solutes of significant size (with a diameter less than 
1.5 nm) able to move as well. The intracellular media con¬ 
tains proteins, particularly the crystallins responsible for the 
high refractive index of the lens. So, the lens is an example 
of a bidomain tissue that has been studied in some detail, 
first in skeletal muscle, then in cardiac muscle and syncytia 
in general. Electrical models of bidomain tissues have been 
developed, and a general approach combining morphology. 
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FIGURE 1 Schematic diagram of lens. (A) The sphere of the lens is 
shown with three landmarks: anterior pole (AP), posterior pole (PP), and 
equator (EQ). ( B ) The control volume in the bidomain model is shown. 
(C) The microstructure of the lens is shown: 1) intracellular region, 2) extra¬ 
cellular region, 3) cell membrane, and 4) gap junction (connections). (D) 
Distribution of the gap junctions between the cell membrane at EQ or AP 
and PP are shown. ( E) A single gap junction that allows the water and 
ion flows is shown. To see this figure in color, go online. 

theory, and experiments has been applied in (1), showing 
how the lens could be studied in this tradition. 

A general approach to bidomain tissues was implemented 

(2) involving detailed measurements of morphology (best 
done with statistical sampling by stereological methods 

(3) ), impedance spectroscopy (4-10) using intracellular 
probes (microelectrodes) that force current to flow across 
membranes to the extracellular baths (1,11-15), electric 
field theory to develop models appropriate to the structure 
(16-20), analyzing the spectroscopic data with the field the¬ 
ory (21,22), and checking that parameters change appropri¬ 
ately (i.e., estimates of membrane capacitance are constant) 
as extracellular solutions are changed in composition and 
concentration (20,23). This work was extended to deal 


with transport by Mathias and co-workers (24-40), and 
computational models of the water flow in the lens were 
later developed in some detail (25,41,42) and exploited 
with great success (reviewed in (43,44); also see (45-49)). 

The original work on electrical models is cited here 
because it provides coherent support, involving a range of 
techniques and approaches, to the general view of syncytial 
tissues used here and in later work. It also shows the range of 
approaches needed to establish a (then) new view of a tissue. 

Mathias (50,51) realized that an asymmetrical electrical 
syncytium would produce convection, in particular in the 
lens (31); he and co-workers systematically investigated 
the flow of water, solutes, and current in the lens, and their 
work is (in our opinion) a model of interdisciplinary 
research, combining theory, simulation, and measurements 
of many types (24-40). Computational models of the water 
flow in the lens were later developed in great detail 
(25,41,42) and compared to the more analytical models. 
These models have been extensively tested, and we are 
fortunate that comprehensive reviews have been written of 
great value to newcomers to the field, particularly (43,44) 
as well as (41,43-49). 

Since the pioneering work on the models of lens 
microcirculation system proposed by Mathias et al. 
(21,51), numerous investigations have been carried out 
(19,20,23,52,53). The microcirculation model has firstly 
relied on a combination of electrical resistance and current 
measurements and theoretical modeling (18,19,54). More 
recently, to provide a better understanding of the electric 
current flow and potential field, the detailed structure of 
lens has been included (41,43-55), describing the asym¬ 
metric biological properties of the lens, and measurements 
of pressure have been made (28,47). Different types of fluid 
flow (56,57) and transport properties of the ions have been 
introduced. Meanwhile, the lens model (55) has been 
extended to simulate age-related changes in lens physiology 
(58) and a variety of physiological processes (26,59-61). 
Reviews of current studies on microcirculation in lens are 
most helpful (26,30,33). Despite this large body of experi¬ 
mental and theoretical work, it is not completely clear 
how they are related to each other. In particular, it is not 
clear how the latest computational models are related to 
the pioneering work and how theoretical analysis is related 
to experimental findings. In this work, we will provide such 
a link. 

Based on the microscale model for semipermeable mem¬ 
branes (62) and the bidomain method (51), we construct a 
mathematical model to ensure that all interactions are 
included and treated consistently. Using asymptotic anal¬ 
ysis, we derive a reduced model, which can be used to obtain 
most physiologically significant quantities except for the 
intracellular pressure. This simplified model can be further 
reduced to the model proposed by Mathias (51) with addi¬ 
tional assumptions that Nernst potentials (which describe 
gradients of chemical potential of each ionic species) and 
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conductance are constant in space. However, we will show 
that neither the Nernst potentials nor the conductance are 
constant. On the contrary, they vary significantly from the 
interior to the surface of the lens. Therefore, both of these 
quantities need to be coupled as part of the solution. 

Our model also shows explicitly that the intracellular 
pressure is decoupled from the rest of the variables. Evolu¬ 
tion has chosen parameters, so the intracellular pressure 
does not affect the other parameters of the lens in a signifi¬ 
cant way. They are robust to variations of intracellular pres¬ 
sure. The evolutionary advantage of this adaptation is not 
clear to us but may be more obvious to other workers with 
a greater knowledge of clinical realities that show how the 
lens becomes diseased (46,48,49,63). Our simplified model 
suggests that all the quantities can be computed without 
knowing the intracellular pressure. On the other hand, we 
need to solve the full model to find the value of the intracel¬ 
lular pressure. Our model is also calibrated by experimental 
data and predicts the effects of gap junctions (28,47) 
described by a “membrane” permeability K in . 

Our new results extend but do not fundamentally change 
previous work on the lens. We strengthen the view that the 
lens provides an osmotic pump to maintain the microcircu¬ 
lation necessary to sustain a living lens for the life of the 
animal. We imagine that similar osmotic pumps create 
microcirculation in other cells and tissues of the body. 

This study is organized as follows. The full model for 
microcirculation of water and ions are proposed based on 
conservation laws in Mathematical Model. In Simplified 
Model, we obtain the leading-order model by identifying 
small parameters in the full model. Based on the boundary 
conditions and partial differential equation analysis, a simpli¬ 
fied version of the leading-order model is proposed and 
compared with the existing models. The model calibration 
and simulation results are shown in Results and Discussion. 
The conclusions and future work are given in Conclusions. 


METHODS 


2. the transmembrane water flux is proportional to the intra/extracellular 
hydrostatic pressure and osmotic pressure differences (i.e., Starling’s 
law, classically applied to capillaries, here applied to membranes 
(64)). In a system like nonideal ionic solutions in which “everything in¬ 
teracts with everything else” (65,66), this statement needs derivation as 
well as assertion. A complete and rigorous derivation can be found in 
(62); 

3. in the rest of this work, the subscript / = in, ex denotes the intra/extra¬ 
cellular space and superscript i = Na + , K + , Cl - denotes the ith species 
ion. 


Then, we obtain the following system for intra- and extracellular veloc¬ 
ities in domain Q — [0, R]: 


1 d 

r 2 dr 


(r 2 M ex u ex ) 


— —M. v L m {P ex — P in + y m ksT(Oi„ — O ex )) 


(la) 


and 


1 d 

r 2 dr 


(/ (“V| ( . v l{ e: 


+ -MinUin)) — 0 , 


(lb) 


where W/ and P/ are the velocity and pressure in the intracellular and extra¬ 
cellular space, respectively, and 0/ is the osmotic pressure with definition 


o„ 


= o = Tc + — 

/ yexi / .in + ,, ! 

v in 


where CJ is the concentration of ith species ion in / space and is the 

density of the permanent negatively charged protein. In this work, we as¬ 
sume the permanent negative charged protein is uniformly distributed 
within intracellular space with a valence of z. Here, J\4/ is the ratio of intra¬ 
cellular area (/ = in) and extracellular area (/ = ex), JA V is the membrane 
area per volume unit, y m is the intracellular membrane reflectance, L m is 
intracellular membrane hydraulic permeability, k B is the Boltzmann con¬ 
stant, and T is temperature. 

As we mentioned before, the intracellular space is a connected space, 
where water can flow from cell to cell through connexin proteins joining 
membranes of neighboring cells, and the extracellular space is narrow, 
with a high tortuosity. The intracellular velocity depends on the gradients 
of hydrostatic pressure and osmotic pressure (41,51,62), and the extracel¬ 
lular velocity is determined by the gradients of hydrostatic pressure and 
electric potential (41,67), 


Mathematical model 

In this section, we present a one-dimensional spherical symmetric nonelec¬ 
troneutral model for microcirculation of the lens with radius R by using the 
bidomain method (17). The model deals with two types of flow: the circu¬ 
lation of water (hydrodynamics) and the circulation of ions (electrody¬ 
namics), generalizing previous bidomain models that deal only with 
electrodynamics. The model is mainly derived from laws of conservation 
of ions and water in the presence of membrane flow between intra- and 
extracellular domains. We note that a similar approach may be useful in 
other tissues with narrow extracellular space, such as the heart, cardiac 
muscle, and the central nervous system, including the cerebral cortex, the 
“brain.” 

Water circulation 

We assume the following: 

1. the loss of intracellular water is only through membranes flowing into 
the extracellular space and vice versa (17); 


K ex d , d 

=- -T c —P ex -k e T c —4 

fi dr dr 


and 


Win 



Y 



(2a) 


(2b) 


where </>/ is the electric potential in the / space, r c is the tortuosity of extra¬ 
cellular region, /x is the viscosity of water, k e is introduced to describe the 
effect of electro-osmotic flow, and /q is the permeability of the intracellular 
region (/ = in) and the extracellular region (/ = ex), respectively. 

Thanks to Eq. 2, Eq. 1 can be treated as an equation of hydraulic pressure. 
Because of the axis symmetry condition, homogeneous Neumann boundary 
conditions are used for pressure at r = 0. At the surface of lens r= R, we set 
the extracellular hydrostatic pressure to be zero, and the intracellular veloc¬ 
ity is consistent with Eq. 2 
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dPex dP itt 
dr dr 

at r 

= 0, 

Pex = 0, 

at r 

= R 

“f {dr Pin ~ 



= L s {P in - y s k B T{O in - O ex )), 

at r 

= R 


where y s is surface membrane reflectance and L s is surface membrane hy¬ 
draulic permeability. 

Ion circulation 

With similar assumptions, the conservation of ion concentration yields the 
following ion flux system: 

=MJ m (4a) 


j‘ = c Uin - iy —c. - o' —c. — d>:„, 

m m ,n dr m 


(5b) 


where D\ is the diffusion coefficient of the ?th species ion in the / space. The 
intracellular medium can be well described by a single diffusion coefficient 
despite it being crowded with solute: current and ions move around the sol¬ 
utes, rather like water flowing around a ship, changing the effective diffu¬ 
sion coefficient and little else, as they do in skeletal muscle (68,69), 
where the cytoplasm (i.e., myoplasm) is even more crowded with ions 
and structures than the lens. 

The Ohm’s law conductance formulation is used to describe the trans¬ 
membrane flux of ions across intracellular membrane and surface 
membrane 


L = (6a) 

j‘s = ^-{<l>in- <t>e*- E ‘), (6b) 


and 


1 d 

r 2 dr 


(r\M ex J' ex + M m J\,)) = 0. 


(4b) 


The ion flux in the intracellular region J\ n and ion flux in the extracellular 
region J l ex are defined as 


d z l p d 

f = C u ex - D l t c — C - D l r c - C —d) py , 

ex ex ex ex c^ e x ex ex d Vex; 


(5a) 


where E 1 — kgT/ez l \og(C l ex /C\ n ) is the Nemst potential (an expression of 
the difference of chemical potential) of the ith species ion. 

In Eq. 6, the intracellular ion conductance g l and surface ion conductance 
G l depends on the ion channel distribution on the membrane (see Fig. 2). 
Based on previous work (17,41,51), we assume that 1) only Na + and Cl - 
can leak between intracellular and extracellular through ion channels inside 
the lens and 2) there is no transmembrane flux for K + between the extracel¬ 
lular and intracellular region, i.e., j= 0. 

Similarly, homogeneous Neumann boundary conditions are used at r = 0. 
At r = R, the extracellular concentrations are fixed, and Robin boundary 


A 



I Na + channel 
K + channel 
I Cl - channel 
Gap junctions 
l Na + /K + pump 



FIGURE 2 (A) Schematic diagram of ion circu¬ 
lation and the distributions of ion channels and 
pumps. The purple line represents the sodium cir¬ 
culation, the light green represents the potassium 
circulation, and the brown line represents chloride 
circulation. The surface epithelial cells ( dark blue 
square ) connect with the intracellular cells {light 
blue hexagon) by the gap junctions {orange rect¬ 
angle). The sodium and chlorine ion channels are 
located on the intracellular membranes, whereas 
the potassium ion channel and sodium-potassium 
ATP pumps are found only on the surface mem¬ 
brane. {B) A schematic diagram of water circula¬ 
tion is shown. Transmembrane water transport is 
through AQPO and AQP1 gap junctions. APQO 
gap junctions are located on the intracellular mem¬ 
branes, and AQP1 is on the surface membrane. To 
see this figure in color, go online. 
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conditions are used for intracellular concentrations because of the trans¬ 
membrane flux and pump: 


4 = 4 = 0, at r = 0 

<4 = c g , 4= ji + a‘, at r = R, 


(7) 


where a 1 is an active ion pump on the surface membrane. Here, we only 
consider the sodium-potassium pump on the surface. The strength of the 
pump depends on the ion’s concentration as in (29,41), 


a Na = 3 /; \ a K = —2—, a a = 0, 


( 8 ) 


where 


C /v 
ft 


4-4,11 c«a +K , \ c *+k k 


c 


~\~Imax2 \ 


£Na \ ^ / qK x ^ 


Cl a + K Na2 ) KCX + Ki 


(9) 


Because of the capacitance of the cell membrane, assumptions of exact 
charge neutrality can easily lead to paradoxes because they oversimplify 
Maxwell’s equations by leaving out altogether the essential role of charge. 
We use the analysis of (70) and thus introduce a linear correction term re¬ 
placing the charge neutrality condition (41,51) without introducing signif¬ 
icant error (see also (71)), 


(i 


v) 



-M v C m ((j) m ~ 4> ex ) (10a) 


and 


V 




M v C m {(j) in - 4> a ), (10b) 


where 7] is the porosity of intracellular region and C m is capacitance per unit 
area. 

Multiplying each ion concentration equation in Eq. 4 with ez, respec¬ 
tively, summing up, and using Eq. 10, the sodium equations are replaced 
by the following equations: 


1 d 

7 dr 



^Pex U ex 





= M v 



4>ex) 



(11a) 


1 d 

7 dr 





= - M 



4>ex) 



(lib) 


with boundary conditions 


' = d<j> in 

dr dr 

4* ex 0 , 



at r = 0, 
at r = R , 


= at r = R, 

< i 

where p in = M,C m /r){<j> in - <j> ex ) + \z\eA in /Vi n and p ex = M,C m /\ - 
V{<t>ex ~ fan) 


8 m 


E*‘- G " = E g ' l t 



In Eq. 11, we define the intracellular conductance a in and extracellular 
conductance o ex as 


= /,;7 C " j 


It is obvious that system 11 might be derived using either Eqs. 4 or 10. 
Therefore, we should drop either Eqs. 4 or 10 when Eq. 11 is used. 

Nondimensionaliza tion 


Because lens circulation is driven by the sodium-potassium pump, it is nat¬ 
ural to choose the characteristic velocity u* n by the pump strength a Na *\ 


a N 

~0*~’ 


where O* — 2(C N 0 a + C*) is characteristic osmotic pressure. Using Eq. lb, 
we obtain the scale of u ex as 


1. * 

— ^0 u ini 


(13) 


where do — A4 ex /Min- With the characteristic values for </>, P, C', chosen as 
ksT/e, fiRu* x /K ex r c , and C N 0 a + C* , we obtain the dimensionless system for 
lens problem as follows (a detailed derivation is given in Appendix C in the 

Supporting Materials and Methods): 

d d 

^ex ~Z~-Pex U 1 , 4*ex 5 (14cl) 

dr dr 

$2 u in = + —Op,, (14b) 

dr dr 




ex — P in) + [Oj n — O ex ), (14c) 


Upx Hit 


(14d) 


^2 Z ' C 'm + = <4(</>,„ - 0J, 


(14e) 
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= ~ S i (&" ~ 0«)> 


(14f) 




l d 
r 2 dr 



<5s 

r 2 

A- 

1 

d 

(r 2 J K ) = 

= o, 

r 2 

dr 

V ex) 


1 

d 

(r 2 J K ) = 

= 0 , 

r 2 

dr 

V in) 

'exUex ~ 


-c 

dr ex 


1 d 

r 2 dr 




(14g) 

(14h) 

(14i) 

(14j) 


dr 


(14k) 


1 d ( 2 ( v—■> / • d ■ d 

J X I ' Pe in fhn U m - /_Pin Z in ~ CT ‘n~r. ( l , in 


<5a d ( i ( v—> ,• ; d ■ d 

= ~7-dr\' I Pe " xfK 'P l " X ~ Y. D e* Z ^fex ~ 


(141) 


and 


e! 



(15d) 


n = 


l P R 

eD* n C *’ 


(15e) 


j\ = Pcqic-D'A-q + z'q-cp, 


(15f) 


Simplified model 

The full model given by system 14 with boundary condition 15 is a coupled 
nonlinear system. In this section, we present a simplified version of the full 
model that captures the main features of the lens circulation. We first obtain 
the leading-order model by identifying the small parameters. Then, by us¬ 
ing boundary conditions and theoretical analysis, the leading-order model is 
further simplified as only one partial differential equation with serial 
algebra equations. 

According to those dimensionless parameters presented in Appendix B in 
the Supporting Materials and Methods, we identify the scale of the param¬ 
eters as follows: 


{<5i ^ <5s} {<5 o 7 ^ 3 } ^0(e 2 ^ 

{<5 2 ,<5 4 A AA} c o(e 2 )- 


If we denote 5g = Df 1 — Df and <5,o = Df 1 — Df a , l = in, ex, it yields 

<5 9 = 0(e 2 ), <5 10 = 0(e). (17) 


with homogeneous Neumann boundary conditions at r — 0 and the 
following boundary conditions at r = 1 : 


Pex = 0 , 

<55 Urn = <5 3 Pin - (0,„ - O cx ) , 

cl = Cl Jf, 

Z K 


[ nvi„ - A), Jfj = o, 


c^cr+cr+i 

<i'„ = 0 , 


P^inPin U in ^ 'P in 1 /h £m a in ,4 


R* 


where 


Z K 


Pin ~ Po + ^6{<t>in ~ Po — \ Z I TT’ 

* in 


Pex ^7 (4*ex 0m) ? 


(15a) 

(15b) 

(15c) 


A priori estimation 

In this section, we provide the priori estimation of the Jfj as follows. By 
using the homogeneous Neumann boundary condition at r = 0, Eq. 141 
yields 






(18) 


From Eq. 18, because Pe in — O(e) and order of <5g, Sg, 610 in Eqs. 16 and 
17, we obtain that 

d d 4> m = 0{e). (19) 


Meanwhile, from Eq. 14b, we can have 

d 


dr 


-0;„ -- ()(e 


( 20 ) 


and in Eq. 14e, we know 


| c " = |(C + C*i)+o(r). (21) 
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With Eqs. 20 and 21 and A in /V, n as constants, we obtain 




( 22 ) 


which means C^ /,G and C^’° are constants and 

dC^° dCf" 


dr 


dr 


, le{in,ex}. 


(30) 


Furthermore, using Eq. 14d and boundary conditions for C c e [ in Eq. 15 
yields 


s-iCl _ /~iNa i s~iK A I ^ I ^in j ^2 

in '“'o ' ^ n 


1+ l-'l ^l + 0 (e z ). (23) 


2 V, 

From the experimental setting of lens (41,51,55), we assume that 

e + cf-!±M£ = o W . 


o o 


2 V- 

^ y in 


(24) 


Therefore, 


In all, we claim that 


Cf: = O(e). 


(25) 


jci - Pe C cl u - D a (—C a 4- z a C a —(h 

J in ~ re in^i n u in U j n l in ' Z ^ in ^ r ^ m 


(26) 


= ^(r). 


By dropping the terms involving these small parameters, the leading or¬ 
der of water circulation system 14a, 14b, 14c, and 14d is as follows: 


u° = --P 0 - 0J-V , 
“ dr ex 1 dr ex 


U n 

= °> 


O u - (f = 0 , 

in ex ’ 


0 0 

U = —U , 
ex in 5 


(27a) 

(27b) 

(27c) 

(27d) 


The leading order of potassium and chloride concentrations satisfy 

1 d 
r 2 dr 


'r 2 J K /) = 0, 


(31a) 


1 d 
r 2 dr 


-ra = o, 


(31b) 


(310 


and 


= <31d| 


where J'f° — PeiC’f — D\(d / drC'f + z'C'fd/drtf) with i — K, Cl and 
/ = in, ex, £ c/ -° = ^log(Cg’°/C^°). 

For the electric potential, using the homogeneous Neumann boundary 
condition at r — 0 and Eqs. 29a, 29b, 30, and 141 yields 


Jr^ - h { PeMn + d9 i C * + 


+— \ PCexpex u ex + ^^fex + ^ W ^f^ex ~ <Tex ^ >ex 


(32) 


At the same time, based on the intracellular equation of potassium 
(Eq. 14j), the homogeneous Neumann boundary condition at r = 0, and 
Eqs. 29a, 29b, and 30, we have 

D ljCl = ( Pe in Ciu in - Df n z K Cf„^^ . (33) 


where the superscript “0” denotes the leading-order approximation. From 
Eq. 27, we deduce 0 Q ex — 0 Q in are constants, and the intracellular and extra¬ 
cellular flow are counterflow. The total charge in the leading-order systems 
is neutral: 


E-^° 


II 

© 

(28a) 

= 0. 

(28b) 


Combining constant osmotic pressure and charge neutrality yields 

O 0 ex (r) = Ol(r) = 2(C N „f(\) + Cf;°(l)), (29a) 


dQl 


dr 


dC^l 

dr 


= 0 , 


Substituting Eq. 32 into Eq. 33 yields 




= ( ( 1 -z K D K A\Pe inUin + -/ D l 
+0(e 2 ), 


( ^8 & ex d(f> e 


dr 


Cl 


(34) 


where we used the fact that p ex = o(e 2 ), d 9 = 0(e 2 ), and dCf /dr — 
— dCl°/dr + 0(e 2 ), le{in, ex}. Because Pe in — 0(e) and <5 8 = 
0(e), in Eq. 34, we claim 


dCi 

dr 


= 0(e). 


(35) 


(29b) Combining Eqs. 32 and 35 yields the leading-order approximation of 

intracellular potential 
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Ui = \pwi - = 0(e), (36) 

(Al U U HI 


where < = 

Similarly, the leading-order approximation of extracellular potential is 


^ ^ ( ll * d r Na,0 , _0 d , o 

T^v r°^ c “ 

= M“(2(0° -0°J - E mfi - E? l » 


(37) 


P° =0, C c, '° = c Wa + c*, c*'° = c*, 

ex ’ ex o ' o i ex o’ 


Pe C K '°u° — D‘ 

re in'-' m 11 in i 


R s 


( d K fi , K r Kfi d 0 \ 
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where i^ 0 = ^log(C^ 0 /C’°)* 

To summarize, the leading-order approximation of system 14a, 14b, 14c, 
14d, 14e, 14f, 14g, 14h, 14i, 14j, 14k, 141, and 15 is given by, in domain 
Q - [0, 1], 


,o d o x d i o 
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In the following, we will further simplify Eqs. 38a, 38b, 38c, 38d, 38e, 
38f, 38g, 38h, 38i, 38j, 38k, 381, and 39 and obtain the relationships be¬ 
tween (jP ex and other leading-order variables by using assumptions concern¬ 
ing the boundary conditions. 

Relation between 0° and </>° ex 


Combining Eqs. 38a, 38d, and 38k and integrating with respect to r yields 
the relation between (jP in and <jP ex as 
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where we used the boundary conditions <j>® x (l) = = 0. 

Relation between P° ex and (jP ex 

By the homogeneous Neumann boundary condition on r — 0 and Eq. 38j, 
we have 
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By Eq. 29b, we can divide Eq. 41 by CfJ’° on both sides, and we get 
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with boundary conditions at r = 1 


(38j) 

(38k) 


Based on the charge neutrality Eq. 28, constant osmotic pressure Eq. 29a, 
and parameters in Appendix B in the Supporting Materials and Methods, we 
define 
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(381) 


Then, combining Eq. 36 and Pe in — O(e), Eq. 42 yields the following 
equation by omitting the higher-order terms: 


Pe„ul . - D c ‘z a 
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Finally, by using the boundary condition, we have the relation between 
extracellular pressure and electric potential as 


pO _ ^ex P e ex&\ Q 

“ ^ Pe a 


(45) 


* “ I „2 / I d r Na,0 I _0 d jfi 

V dr 1 1 ' V C “ + ex dJ ex 


(51a) 
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Expression of E Na 

Based on potassium equation and relation in Eqs. 40 and 45, we have the 
expression for Cf n and C* as 
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Based on Eq. 28, we can get 
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Extracellular electric potential system 

By Eqs. 40 and 45, we have 0 in as 


\ °in re ex &i n J 


(49) 


The value 0^(1) is determined by the boundary condition of 0? in Eq. 
39, where 

-M* [ (2(0° -0°,) - E Na '° - E clfl )s 2 ds = a"*, (50) 
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where we use 
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To summarize, we obtained the simplified model of system 38a, 38b, 38c, 
38d, 38e, 38f, 38g, 38h, 38i, 38j, 38k, 381, and 39 as follows: 
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(52) 


Under the same assumptions in (51)—for example, uniform diffusion 
constants for all ions and constant Nemst potential—our simplified model 
system 5 1 recovers the model proposed by Mathias. The main contribution 
here is that we remove the assumptions that Nemst potentials and effective 
conductance should be constants. By using the relationships between ion 
concentrations and external potential, we obtain the space-dependent 
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Nernst potential, which yields a much better approximation to the full 
model (see Fig. 4). 

RESULTS AND DISCUSSION 

In this section, we present numerical simulations using both 
the full and simplified models. The finite volume method 
(70) is used to preserve mass conservation of ions. The 
convex iteration (72) is employed to solve the nonlinear 
coupled system. The numerical algorithm is implemented 
in MATLAB (The MathWorks, Natick, MA). 

Model calibration: Membrane conductance 
effects intracellular hydrostatic pressure 

In this section, we first calibrate the full model by compar¬ 
ison with the experimental data to study the effect of con- 
nexin on intracellular hydrostatic pressure. 

Intracellular hydrostatic pressure is an important physio¬ 
logical quantity (63). In some studies (28,47), the authors 
showed that the connexin (gap-junction; see Fig. 1, D and 
E) conductance plays an important role in the microcircula¬ 
tion of lens. It is said that if the intracellular conductance 
Kj n /Hi n in lenses is approximately doubled, the hydrostatic 
pressure gradient in the lenses should become approxi¬ 
mately half of the original one. In this section, we calibrate 
our model. We choose a value of the intracellular conduc¬ 
tance (k,■„//!„,) that correctly calculates the experimental re¬ 
sults in (28,47). 



In Fig. 3 A, the value k” = 4.6830 x l(T 20 /m 2 (black 
line ) yields a good approximation to experimental data 
0 blue markers). When the conductivity of the connexins is 
doubled so that parameter value K in is 2k'/ (in the lens of 
mice Cx46 knock-in [KI] lens) as in the experiments 
(28,47), which doubled the conductivity of the connexins 
by using Cx46 KI mice lens, our model (black dot) can 
also match the experimental data (red markers): the intracel¬ 
lular hydrostatic pressure drops to half. This result shows 
that our full model can correctly predict the effect of perme¬ 
ability of membrane on hydrostatic pressure. 

Interestingly, Fig. 3, B—D show that other intracellular 
quantities and extracellular ones (Supporting Materials 
and Methods) are insensitive to increases in the permeability 
by a factor of 20, even to 2(k" n . The reason for this can be 
explained by using our simplified system 51 . If the variation 
of intracellular conductance still keeps the <52 to be a small 
quantity in the dimensionless system 14, our simplified 
model will be still valid. In the simplified model, all the 
quantities except intracellular hydrostatic pressure are 
related to the extracellular electric potential. However, the 
extracellular electric potential will not be affected by the 
change of the intracellular conductance because Eq. 51a 
does not involve intracellular conductance. 

Full model versus simplified model 

In this section, we compare the full model 14a, 14b, 14c, 
14d, 14e, 14f, 14g, 14h, 14i, 14j, 14k, 141, and 15 with the 




FIGURE 3 Comparison between different K in . (A) Comparison between simulation results and experimental results. The experimental data of dog, rabbit, 
and rat come from (47). Mice and Cx46 KI mice come from (28). According to (28), the Cx46 KI mice lens has twice the number density of lens gap-junction 
channels compared to mice. The parameter K™ n — 4.6830 x 10 -2O /m 2 , and radius is written in dimensionless units for different species. ( B ) Space distribution 
of intracellular osmotic pressure (O in ). (C) Space distribution of intracellular electric potential (</>/„). ( D ) Space distribution of intracellular water velocity 
(u in ). To see this figure in color, go online. 
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simplified model 51 and Mathias model in (51). The nu¬ 
merical results of full model (black lines ) in Fig. 4, A-C 
suggest that the variations in intracellular electric poten¬ 
tial, extracellular conductance, and Nernst potential of 
Cl - are rather small. The assumption of constant 
values for those variables (potential, extracellular conduc¬ 
tance, and Nernst, i.e., chemical potential of Cl - ) in the 
Mathias model (shown as red dash-dot lines ) is reason¬ 


able. However, the Nernst potentials of sodium and 
potassium (Fig. 4, D and E ) have large variations because 
of the effect of sodium-potassium pump. Our simplified 
model (blue dash lines ) describes these variations with 
small errors. The comparisons for extracellular pressure, 
velocity, and potential (Fig. 4, F—H) confirm that our 
simplified model yields good approximations to the full 
model. 



FIGURE 4 Comparison of electroneutral and simplified and the Mathias model in (51). (A) Space distribution of intracellular electric potential (</>,„). ( B ) 
Space distribution of extracellular conductance {(T ex ). (C) Space distribution of chloride nernst potential ( E cl ). ( D) Space distribution of sodium nernst 
potential (E Na ). (E) Space distribution of potassium nernst potential ( E K ). ( F) Space distribution of extracellular pressure (P ex ). (G) Space distribution of 
extracellular velocity ( u ex ). ( H) Space distribution of extracellular electric potential ((p ex ). To see this figure in color, go online. 
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CONCLUSIONS 

In this work, we propose a bidomain model to study the 
microcirculation of lens. We include a capacitor in the rep¬ 
resentation of the membrane, and so our model is consistent 
with classical electrodynamics. Consistency produces a 
linear correction term in the classical charge neutrality 
equation. This full model is calibrated by comparison with 
the experiment studying effect of connexin on hydrostatic 
pressure. By changing only the intracellular membrane 
conductance (strength of connection), our model could 
match the two experimental results with different connexins 
very well. Our model is capable of making predictions about 
the circulation of the lens. Furthermore, the numerical sim¬ 
ulations show that the velocity, potential, and osmotic pres¬ 
sure in the intra- and extracellular spaces are not sensitive to 
increasing intracellular membrane conductance. 

Based on the asymptotic analysis, we proposed a simpli¬ 
fied model that allows us to obtain a deep understanding of 
the physical process without making unrealistic assump¬ 
tions. Our results showed that the simplified model is a 
good approximation of the full model in which Nernst po¬ 
tentials and conductivity vary significantly inside the lens. 

Our model allows calculation of variables that determine 
the role and life of the lens as an organ. Particularly impor¬ 
tant are the factors that determine the transparency of the 
lens because that is the main function of the organ. The 
dependence of the size of the extracellular space, and thus 
the pressure in extracellular and intracellular spaces and 
the difference between those two, is likely to be an impor¬ 
tant determinant of transparency. One imagines that 
swelling of the extracellular space will scatter light, partic¬ 
ularly because the swelling is likely to be irregular (in a way 
our model does not yet capture). Changes in the osmolarity 
(i.e., activity of water estimated by the total concentration of 
solutes) are likely to be important as well. 

This hydrodynamic bidomain model can point the way to 
dealing with other cells, tissues, and organs in which current 
flow, water flow, and cell-volume changes are important. 
These include the kidney, the central nervous system (where 
the narrow extracellular space poses many of the biological 
problems facing the lens), the t-tubular system of skeletal 
and much cardiac muscle, and so on. We show that a mathe¬ 
matically well-defined model can deal with the reality of bio¬ 
logical structure and its complex distribution of channels, etc. 

Conservation laws applied to simplified structures are 
enough to provide quite useful results because they were 
in three-dimensional electrical problems of cells of various 
geometries (16) and syncytia (3,23). The exact results are 
analyzed with perturbation methods, described in general 
in (73), and these methods allow dramatic simplifications 
without introducing large or even significant errors. It is 
as if evolution chose systems in which parameters and struc¬ 
tures allow simple results and in which parameters can con¬ 
trol biological function robustly. 


Of course, we only point the way. Additional compart¬ 
ments and additional structural complexity will surely be 
needed to deal with the workings of evolution. But these 
can be handled in a mathematically defined way, yielding 
approximate results with clear physical and biological inter¬ 
pretation. By combining the multidomain model and mem¬ 
brane potential dependent conductances, one can model 
depolarization induced by extra potassium in the lens 
(53,55) and the cortical spreading depression problem 
(71,74,75). The ultimate goals will be 1) to provide as 
much precision in the mathematics and physics as we can, 
starting from first principles (62), and 2) to provide a general 
basis for treatments of convection in other tissues that 
involve microcirculation. Computational models of these 
are not in hand and may be hard to construct because so little 
is known of those systems compared to the lens. With what 
we have learned here, we hope a general mathematical 
approach and model of the type we present here may be con¬ 
structed and prove helpful in other systems with narrow 
extracellular spaces that are likely to need microcirculation 
to augment diffusion, such as cardiac and skeletal muscle, 
kidney, liver, epithelia, and the extracellular space of the 
brain. 

SUPPORTING MATERIAL 

Supporting Materials and Methods, one figure, and two tables are avail¬ 
able at http://www.biophysj.org/biophysj/supplemental/S0006-3495(19) 
30134-1. 
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Appendix A. Model Parameters 


Parameters 

Mathias [1] 

Malcolm [2] 

Parameters 

Mathias [1] 

Malcolm [2] 

R 

1.6 x 10~ 3 m 

1.6 X l(r 3 m 

Lm 

3.75 x 10~ 13 m/ [Pa ■ s) 

1.34 x 10- 13 m/(Pa-s) 

Ain / Vi„ 

78 mM 

78 mM 

L s 

3.75 x W- 13 m/(Pa s) 

8.89 x 10 - 13 m/(Pa s) 

fNa 

107 mM 

107 mM 

Mi n 

0.988 

0.99 

c K 

'~o 

3 mM 

3 mM 

A'iex 

0.012 

0.01 

Cm 

- 

1 x 1(T 2 F/m 2 

M„ 

6 x 10 5 /m 

5 x 10 5 /m 

D Na 

‘-'ex 

- 

1.39 x 1CT 9 m 2 /s 

T 

- 

310 K 

d k 

ex 

- 

2.04 x 10~ 9 m 2 /s 

K 

1.72 x 10 - s m 2 /(V-s) 

1.45 x 10~ 8 m 2 /(V-s) 

D a 

‘-'ex 

- 

2.12 x 10~ 9 m 2 /s 

kB 

1.38 x lO^ 23 J/K 

1.38 x 10- 23 J/K 

D Na 

- 

1.39 x 10“ n m 2 /s 

Kn 

- 

1.6154 mM 

D f„ 

- 

2.04 x 10- 11 m 2 /s 

Kk2 

- 

0.1657 mM 

Dg 

- 

2.12 x 10 -11 m 2 /s 

K-Nal,Na2 

- 

2.3393 mM 

e 

1.6 x 1(T 19 A ■ s 

1.6 x 10~ 19 A ■ s 

>7 

0.988 

0.99 

g N ‘ 

2.2 x 1(T 3 S/m 2 

2.2 x 10~ 3 S/m 2 

Kex 

1.141 x 10~ 16 m 2 

1.33 x 10~ 16 m 2 

g a 

2.2 X 1(T 3 S/m 2 

2.2 x 10~ 3 S/m 2 

Kin 

- 

9.366 x 10~ 19 m 2 

G K 

2.1 S/m 2 

2.1 S/m 2 

7m,s 

1 

1 

b 

2.3 x 1( r 2 A/m 2 

- 

%c 

0.16 

0.16 

hnaxl 

- 

0.478 A/m 2 

¥ 

7 x 10~ 4 Pa ■ s 

7 x 10~ 4 Pa ■ s 

hnax2 

- 

0.065 A/m 2 

z 

-1.5 

-1.5 
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Appendix B. Dimensionless Parameters and Scales 


The following dimensionless parameters' value and scales calculation based on values in [2] 


Scales / Parameters 

Value 

Parameters 

Value 

a N «* 

6.9 X 10 - 2 A/m 7 

*0 

l 

99 

C* 

110 mM 

Si 

1.2031 x 10- 1 

O* 

220 mM 

h 

6.861 x 10~ 3 

P * 

16.937 KPa 

S3 

2.9894 x 10~ 2 

K, 

3.2506 nm/s 

Si 

3.5323 x 10- 5 

<x 

3.2181 pm/s 

Ss 

4.3022 x 10~ 3 

<P * 

26.7 mV 

Sf, 

1.2745 x 10~ 5 

D ex 

3.392 x 10 10 m 2 /s 

s 7 

1.2617 x 10~ 3 

DU 

2.12 x 10~ n m 2 /s 

s s 

1.6162 x lO^ 1 

Pe e x 

1.5180 

Si 

12.5 

110 

Pern 

2.4533 x 10- 1 

<$10 

3.443 x 10- 1 

jjNa 
in,ex 

0.6557 

Si 1 

3.77 x 10~ 2 

in,ex 

0.9623 

Po 

117 

110 

m 

in,ex 

1 

Mf 

3.3859 x 10- 1 

R s 

4.00 x lO^ 1 

M? 

2.095 


2 








Appendix C. Non-dimensionalization 

In this section, we derive the dimensionless model based on the lens, which has been widely stud¬ 
ied. The major ions we considering here are sodium ( Na + ), potassium ( K + ) and chloride (C/ _ ) and the 
sodium-potassium pump which distributed on the surface of the lens. Although we restrict ourselves in 
this particular problem, the following procedure can be applied in a wide range of practical problems in 
biological syncytia. 

Appendix C.l. Water circulation 

In the following, we assume the typical length of lens is R. The fluid system is driven by the osmotic 
gradient, which is generated by the sodium-potassium pump on the surface. In Eq. 7, the strength of 
sodium-potassium pump at surface depends on the ion's concentration, which leads 

a Na —3—, a K = - 2^, a cl = 0, (C.l) 

e e 


where 


Ip I max] 


c 


■Na 


c!i a + K 


Na 1 


c* 


cf- 


Kki 


Imax2 


c 


Na 


t-Na i v 

+ K Nn2 


C* 

v —n 


cf 


Kk2 


(C.2) 


We assume that the velocity at surface determines the characteristic velocity scale for the problem. We have 
ion fluxes in the intracellular, extracellular region in Eq. 5 and trans-membrane source of ion in Eq. 6 for 
ion Na + ,K + , Cl~. 

At boundary of the intracellular space, due to the ion pump in Eq. C.l and assumption of conductance at 
surface that G Na — G cl — 0 [1, 3], we have 


,-Na _ Na jK 
‘in u ' Jin 


]s+a\ = 0 . 


(C.3) 


Since g K — 0 inside of the lens, we obtain 

jf + n K — 0. (C.4) 

This assumption obviously will have to be replaced in applications to other tissues, with a less particular 
distribution of channel proteins. 

By the conservation of fluxes for each ion in Eq. 4, we get 


Ji„ — ^oJex' i = Na,K,Cl, 


(C.5) 


3 






where <5 q = Therefore, Eq. C.3 becomes 


-6 0 J? x “ = a N “, -SojX= 0, -S 0 jg=0. (C.6) 

Adding up all three fluxes in Eq. C.6 and since in the extracellular region each ion diffusion coefficient are 
at the same level of approximation, i.e. 

Dg X — O (D ex ), i — Na, K, Cl, (C.7) 


and based on Eq. 10 , we get 


i r n ^ r\ i ^ex^c^O d Na 

Oex^in T OoD ex T c -j-O ex H ; V Pe x ~7~(pex — 0 
dr kgT dr 


(C.8) 


The strength of the ion pump a Na depends on the ion concentration in Eq. C.2 . We choose the scale of a 


Na 


is a Nn * based on an experimental estimation [1], Using Eq. C.8, we take the scale for O ln ex and u m to be O* 


and u* n as 


O* =2(C 0 Na + C5V u* n = a 


Na * 


in q * 

By mass conservation expressed in Eq. 1, we naturally get the scale of u ex as 


(C.9) 


* f—1 * 

u ex ~ ^0 u in- 


(C.10) 


Furthermore, <p * = is used for the scale of electric potential (pj n and cf) ex . For the extracellular velocity 
in Eq. 2, we have 

k _ d ~ hr d ~ 

(C.U) 


k B T d 


u ex u ex — „,T> r cPex jyPex k e T c l Jf'k ex ' 


}iR 

We think the ^ P ex term balance the velocity u ex . The scale for extracellular pressure P* x is then choose 


P* — 


P Ru *ex 

KexTc 


Therefore, we get 


d ~ d ~ 

U ex — ~F^Pex ~x^(pe X / 

dr 


dr 


(C.12) 
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where T . For the intracellular velocity, we have 


K inPj n d ~ Kj n 'Y m kBT(y r d 

/<R dr + /.(R (lr ' r 


(C.13) 


We claim term dr p m and balance at the same level. Therefore, we choose the same scale for the 

intracellular and extracellular pressure, namely. 


p* — p* — p* 

1 1 in L ex' 


ThenEq. C.13 becomes 


d ~ d ~ 

^ 2 11 in ~ dfOinr 


(C.14) 


where 


^2 — 


}iRu* 


Kin7mk B TO 


:> S 3 = 


j m k B TO* 


In all, the fluid system Eq. 1 becomes 


i’Y — U in 


1 d 


5 ^a?\ r 


— ^3 (Pex Pin'j + ^ 0 Bl O ex ^j 


(C.15) 


with boundary condition 


where 


Pex — 0 , 

^5Min Win (^in Oex ) / 






7/ <% = 


RM v L m j m k B TO* L s j s k B TO* 

Appendix C.2. Ions circulation 

The velocity scales and diffusion coefficients in the extracellular and intracellular space are at different 
levels of approximation in our approach. In the following, we put the characteristic diffusion coefficients 
at intracellular and extracellular region and scale of concentration as 


D* t~\CI —- pP pjC/ p* pNa i r-'K 

u ex — u ex T cr — u [ n , ^ ~ ■ 


5 









In this way, we get Peclet number in the extracellular and intracellular and dimensionless Nernst potential 


as 


p e . — U in R p e _ U ex R gz — — log | — 

£w “ Df ' Pe « - D* ' E - l0§ CL , 


Because g Na — 0 inside of lens, we have K + system as in Mathias's model [1], 


— — (r 2 f Pe C K u — D K f —C K 4- z K C K 

r 2 dr \ l Uex l dr ex + Z ^ ex dd e 

— —(i 2 - (p e - C K il • - D k ( — C K + - K C K — 

^2 l ' l 1 t in K -'in U in l r—^in \ 6 '—in r— 


dr 


dr 7 


-0, 

- 0 , 


(C.16) 


with boundary condition 


pK — pK 
^ex '~o > 


Pe- c K d -n K I —C K + z K &—d>- 

1 Zin^in u in \ flyin'* ^in 




and Cl system as 


Ji ( r2 (PeexC%u ex - Dg (jf?x 


7 o pa “ X 
Z Lex dr x<!X 


1 d 
r 2 dr 


r2 I p pa r _ pci ( ArCl, -Cl pci d 

1 t in'-in U in l ^ 


dr 


dr* in 


1 d 


^r 2 dr 1 T 


M e v x 

vCl 


2 1 Pe C c, u - D C1 ( —C a +z cl C a d 

1 tex^ex u ex u e x l 


dr 


<pex £ 


Cl 


(C.17) 


with boundary condition 


cfx — Co 1 " + Of + $7 ((pin — (pex) / 


Pe- C°u- - D° [ — C a + - a C a —d- 1 - 0 


where 


_ G K k B TR K _ a K R ~ _ M v g a k B TR 2 
s e 2 D* in C* ’ D * n C*' v M ex e 2 D* x C* ’ 


■Mex^tx 
KA. n* ' 

jv hn u in 


The concentration of Na + can be solved from the following equations 


Y2 Z ‘ C in dZirf-^Se (t pi n - (pex) , 
* i 


E Z ' Ct ex = -<$7 (<Pm ~ (Pex) , 


(C.18) 
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where 


(C.19) 


r _ MvCmksT „_ MvCmkftT 

6 _ e 2 C*?; ' 7 ~ e 2 C*(l — rj) ’ 

From Eq. 11 and use the fact z ,Vn = z K — 1 and assumption that g Na — g CI and G Na — G cl — 0, we 


have 


,• d 


r ( Pe ex p ex Uex - Y^D'exZ 1 jf Cl ex ~ VeXj~<pex ) ) = M e f [2 ( p in - <p ex ) - E Na - E CI J , 


‘dr 


1 d 
T 2 dr 

\ \ L 

Wd¥\ r \ ^ ein P inUin ~ E Din z tf C in ~ a ‘n j-Qii 


(C.20) 


1 d 


d 


~ ^ 8 ?2 dr ( r ( P e exPex u ex Yl^ex 2 '^-^ex a ex j^P* 


d? 


df 


17 ex I I / 


with boundary condition 


f pex — 0/ 

\ P e inPin u in ~ Yh^in z ^~Qn ~ ^in^-Pin ~ ('Pin ~ E j ~E Ip, 


(C.21) 


where 


and 


Pin — |z| + ^6 {pin pex) / Pc* — <^7 (pex Pin ) / Ip — p* / 

''in 


= E D /n ( zI ) 2C in^ ^ = E Dk'(z ! ) 2 C ( ' 
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Figure D.l: Comparison between different K;„. 
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